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ABSTRACT 

We propose that there exists a class of transient sources, "squeezars" , which are stars caught in highly 
eccentric orbits around a massive (m < 10 8 M Q ) black hole (MBH), whose atypically high luminosity 
(up to a significant fraction of their Eddington luminosity) is powered by tidal interactions with the 
MBH. Their existence follows from the presence of a mass sink, the MBH, in the galactic center, which 
drives a flow of stars into nearly radial orbits to replace those it has destroyed. We consider two limits 
for the stellar response to tidal heating: surface heating with radiative cooling ("hot squeezars") and 
bulk heating with adiabatic expansion ("cold squeezars"), and calculate the evolution of the squeezar 
orbit, size, luminosity and effective temperature. The squeezar formation rate is only ^0.05 that of tidal 
disruption flares, but squeezar lifetimes are many orders of magnitude longer, and so future observations 
of squeezars in nearby galaxies can probe the tidal process that feeds MBHs and the effects of extreme 
tides on stars. The mean number of squeezars orbiting the Galactic MBH is estimated at 0.1-1. 

Subject headings: black hole physics — galaxies: nuclei — stars: kinematics 



1. INTRODUCTION 

A massive galactic black hole is a mass sink that drives 
an inflow of stars from its radius of influence, r^. When 
the MBH mass m is small enough (m < 10 s M Q for solar 
type stars), the tidal disruption radius r t lies outside the 
event horizon and the star is disrupted before falling in. 
Prompt disruption occurs when stars are scattered into 
highly eccentric (loss-cone) orbits with periapse r p < r t . 
The accretion of the debris may be observed as a short- 
lived (~lyr) luminous "tidal flare" (Frank & Rees 1976). 

Stars with r p > r t narrowly escape after undergoing ex- 
treme tidal distortion, mass-loss, spin-up and subsequent 
mixing, which may affect their evolution and appearance. 
The near-miss event rate is comparable to that of tidal 
disruptions (Alexander & Livio 2001). Most of these stars 
avoid tidal capture by being deflected to wider orbits, or 
by missing the MBH due to its Brownian motion. Such 
"tidally scattered" stars eventually amount to a few per- 
cent of the stellar population within ru- The rest are 
tidally captured and gradually spiral in as the tides con- 
vert a small fraction of their orbital energy into heat each 
peri-passage (Alexander & Hopman 2003, AH03). 

The orbital energy a star must lose to circularize far 
exceeds its own binding energy. A "squeezar" , a tidally 
heated star orbiting a MBH, is ultimately disrupted by 
expanding beyond its Roche lobe or by radiating above 
its Eddington luminosity (Rees 1988; Novikov, Petchik 
& Polnarev 1992). Squeezars directly trace the tidal dis- 
ruption process, which is important for feeding low-mass 
MBHs. The squeezar phase lasts orders of magnitude 
longer than a tidal flare, and so observations of squeezars 
in nearby galaxies, in particular IR observations of the 
Galactic Center (GC), can probe the tidal disruption pro- 
cess and the effects of extreme tides on stars. 

2. SQUEEZAR EVOLUTION 



The evolution of a squeezar reflects its structure and the 
way its mechanical and thermal properties respond to tidal 
heating. We approach the challenging problem of model- 
ing the evolution by considering two simple scenarios that 
likely bracket the range of possible responses: (1) Surface 
heating and radiative cooling ("hot squeezar", HS), where 
the tidal oscillations dissipate in a very thin surface layer 
that expands moderately and radiates at a significantly in- 
creased effective temperature, T* (McMillan, McDermott 
& Taam 1987). (2) Bulk heating and adiabatic expansion 
("cold squeezar", CS), where the oscillations dissipate in 
the stellar bulk and cause a large quasi adiabatic, self- 
similar expansion at constant T* (Podsiadlowski 1996). 

We consider the orbital decay of a star of mass M* <C m 
and initial radius i?* that is tidally captured by a MBH 
in a region where the MBH dominates the potential (Ke- 
plerian approximation). We denote by a tilde quanti- 
ties in dimensionless units of G=M i ,=R- k =\. Thus time 
is measured in terms of the stellar dynamical time = 
y/Rl/GM*™ 1600s (i?*/i? ) 3 / 2 (M*/M )- 1 / 2 and energy 
in terms of the stellar binding energy (up to a factor), 
E* = GMt/R* ~ 3.8 x 10 48 erg {M*/Mq) 2 {RJ Rq)- 1 . 

The tidal disruption radius, rj ~i?(f)m 1 / 3 (i? is the stel- 
lar radius) , lies outside the event horizon of a non- rotating 
MBH when m < (c/ v 7 ^) 3 (c is the speed of light) . The semi- 
major axis, period and eccentricity of a Keplerian orbit are 
related to the orbital energy E by 

a= -m/2E, P = 2ir y / a 3 / (1 + m) , e = l-r p /a. (1) 

If orbital angular momentum is conserved, the circu- 
larization radius is r c = 2r p . Conversely, if the angular 
momentum transfered from the^orbit to the tides follows 
the impulsive relation A J = AE/Q P , where Q p is the or- 
bital angular velocity at periapse (e.g. Kumar & Quataert 
1998), then r p = r c = const. In either case r p ~ const in the 
early stages of circularization. We assume here r p = const, 
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so the orbital energy the star has to lose in order to circu- 
larize from an E^O orbit is E c — m /2r p = fh 2 / 3 /2b > 1. 

The tidal energy extracted from the orbit in a single 
peri-passage is given to leading order in the linear multi- 
pole expansion by (e.g. Press & Teukolsky 1977) 



AE t = T 2 ( V ) R 5 /b 6 



(2) 



where b = r p /r t (t — 0) and the R 5 term accounts for pos- 
sible stellar expansion. The 2nd order tidal coupling co- 
efficient Ti depends on the stellar structure, on e, and on 
the dimensionless transit time, r\ = f\j 2 / ' \J\ + m ~ b 3 ! 2 
(rh^> 1). The linear mode analysis is formally valid for 
rj> 1, which holds here. We assume here that e = 1 since 
the P-e relation, 1 — e ~ (2tt / P) 2 / 3 b, indicates that e~ 1 
down to F~ 0.1 yr (e.g. at that point 1 — e = 0.02& for 
a solar type star), and most squeezars are disrupted well 
before reaching such short periods (§3). 

The tidal coupling changes as the stellar structure and 
spin evolve. The analysis below is greatly simplified by 
assuming that T 2 is constant in time. This is suggested 
by the fact that the squeezar phase coincides with the ini- 
tial stages of circularization and synchronization. In the 
case of bulk heating (§2.2), the star expands roughly self- 
similarly and so the structural changes can be accounted 
for by the R 5 term, while in the case of surface heating the 
expanding layer involves only a minute fraction of the stel- 
lar mass and the stellar bulk is unaffected (§2.1). We do 
not consider possible resonances between the orbital pe- 
riod and the tidal oscillations (Novikov et al. 1992). This 
is not relevant for CSs, where the tidal energy is carried by 
high order modes and quickly dissipated. For HSs, we note 
that small orbital perturbations (5P, 5E) by other stars 
will randomize the phase between the orbit and the tidal 
oscillations as long as t/P <SP/P=35E/2E ~3P/2t r , 
where r <t+ is the width of the resonance, and t r is the 
2-body relaxation time. E.g., for r = O.lt* <~ 100s and 
t r = 10 9 yr, resonances are suppressed for P > 50 yr. Fur- 
ther analysis is needed to validate these assumptions. 

Once the details of the tidal energy deposition and the 
stellar response are specified, Eqs. (1) and (2) can be 
evolved numerically orbit by orbit. We now derive approx- 
imate analytic solutions for the evolution of squeezars. 

2.1. Hot squeezars: surface heating and radiative cooling 

The place where the tidal oscillations thermalize de- 
pends on the modes that carry the energy. In stars with 
large convective envelopes, it is mainly carried by / and p 
modes that, absent non-linear couplings to higher modes, 
dissipate in the outermost layers (e.g. McMillan et al. 
1987). We adopt this as a limiting case, and assume that 
R = 1 in Eq. (2) since the expansion involves only a thin 
surface layer with a minute fraction of the stellar mass. 
The implied assumptions are that the oscillations occur in 
the unchanged stellar bulk just below the surface layer in 
which they dissipate, and that the oscillatory modes are 
unaffected by the modified surface boundary conditions. 

The orbhaj_evolutiori_of a^HS^is then derived^ from the 
relation dP/P~(dP/dE)AE t /P (assuming \ AE t /E\<£l, 
as is the case here), whose solution is 



where P$ and ao are the initial values and where 

t = [4^ 2 P /(1 + m)] 1/3 m/A£ t = 2E Q /(dE/dt) Q (4) 

The formal result P(t n ) = is an artifact of the assump- 
tion AE t ~ const. If the star could avoid expansion and 
tidal disruption, then ultimately AE 1 — >0 (Hut 1980) and 
P ^2-Kb 3 / 2 . Because the final stages of the stellar expan- 
sion are rapid, t also estimates the total time to disruption 
for b> 1.5. The number of orbits completed by time t is 

N olh = Jo dt'/P(t') = (t /2P ) [(1 - tfto)- 2 - 1] . 

The energy in the stellar oscillations is dissipated and 
radiated over a timescale t^, which can span several orders 
of magnitude depending on the stellar structure and the 
nature of the excited mode, but usually r d > P. Here 
we adopt a typical value of r d ~ 10 yr (Ray, Kcmb- 
havi & Antia 1987; McMillan et al. 1987). We assume 
that the tidal luminosity L t (the luminosity in excess of 
the initial luminosity L+) that is released by one peri- 
passage decays exponentially, L t = (AE t /Td) exp(— t/rd), 
where t is the time after periapse. The luminosity from 
many successive peri-passages is given by the relation 
L t (t + P) = [L t (t) + AE t /T d ] exp(-P/?a). This can be 
solved approximately by substituting P = const, since most 
of the luminosity comes from the recent peri-passages that 
have roughly the same period. The luminosity is then 

Lt = (AE t /r d ) [l-expt-i/r.,)] /[exp(P> d )-l] , (5) 

for the initial condition L t (0) = 0. At later times, when 
t^>Td and P <C Td, L t — > AE t /P irrespective of r d . 

Numeric models of tidal heat dissipation in a thin sur- 
face layer (McMillan et al. 1987; Podsiadlowski 1996) in- 
dicate that the stellar interior is essentially unaffected and 
that in the limit L t 3> 1 (L t = L t /L±), the binding energy 
required to raise the expanded layer above the stellar sur- 
face can be empirically expressed as a power-law of the 
tidal luminosity, AE± cxl — 1/Roc Lf. Therefore 



R- 1 = 1 - (L t /L ) a {L\ » 1) , 



(6) 



where Lq is the terminal luminosity for which R^oo. The 
results of McMillan et al. (1987) indicate that a~0.2 for a 
0.8 Mq main sequence (MS) star, 0.4 for a 1.5 M star, 

and that 2x 10 3 <L < 10 4 , depending on the depth of the 
heated layer (smaller values correspond to deeper layers). 

The terminal luminosity Lq is of the order of the stellar 
Eddington limit L E = 3.2 x 10 4 (M^/M Q )( J L^/i Q )- 1 . The 
stellar expansion is truncated when R = b, and the tidal 
disruption radius overtakes the star. From that point on, 
the stellar lifetime is limited by mass loss at periapse. Here 
we conservatively assume prompt destruction when R = b. 
The change in the effective temperature, T = T/T+, is 



(1 + L t )/R 2 



(7) 



P = P (l-t/t ) 3 



a = a (l - t/t 



(3) 



Since L t < Lq while R diverges, there exists a maximal 
temperature, which is attained when P=l+l/2a, or equiv- 
alently when L t /L = (1 + 2a)- 1 /" (L t > 1), 

maxf ~ V2^(l + 2a)-^ +2a ^ ia ty i ~0.4lJ /4 . (8) 
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For the a and Lq range considered here, 2.3 < max T< 4.6. 
For T+ > 10 4 K, the change in the stellar K-band magnitude 
can be roughly estimated by the Rayleigh-Jeans limit, 

AK = K - K+ ~ -2.51og 10 ( J R 2 f) , (9) 

where K* is the star's initial magnitude. The actual mag- 
nitude is ^0.5 m brighter for a solar-type star (by numeric 
integration of the blackbody spectrum; the values quoted 
below are exact). The X-luminosity increases with the 
stellar radius and so is largest just before disruption. At 
maximal temperature, — 4.1 m < AK <— 3.0 m . 

These analytic expressions fully describe the evolution 
of a HS given initial orbital parameters (6 and Po) and 
stellar type (M*, R+, L±, T+ and K+). T 2 can be evaluated 
numerically for a given stellar model. Here we represent a 
typical star by the Sun, and calculate AE t for a detailed 
solar model (Alexander & Kumar 2001). Figure 1 confirms 
that the exact orbit-by-orbit evolution of HSs (§2) is well 
approximated by Eqs. (3-9). 

2.2. Cold squeezars: bulk heating and adiabatic expansion 

Stars subjected to high amplitude tidal deformations 
may dissipate the tidal energy in their bulk via non-linear 
mode couplings, which redistribute the energy among an 
infinitude of high order modes (Kumar & Goodman 1996). 
These dissipate very quickly, up to 10 6 times faster than 
low order modes, and so bulk heating is effectively instan- 
taneous. Numeric modeling of bulk heated stars (Podsi- 
adlowski 1996) indicates that they expand quasi adiabat- 
ically and self similarly on a dynamical timescale while 
maintaining their original effective temperature. 

The tidal heat is stored in the stellar binding energy, 

AE* = (8(1 - l/R) , (10) 

where the dimensionless factor j3 depends on the stellar 
mass concentration and heat capacity ratio. Typically 
(3 = 0.75 for a MS star (e.g. Cox & Giuli 1968). The 
tidal heat is radiated from the expanded surface, and so 
the luminosity peaks just before tidal disruption, when 
L = b 2 . Since typically L t < 10, CSs are far less luminous 
than hot ones, but because AK is dominated by the stellar 
expansion, they brighten by up to — 5 log 6 ~ — 2.5 m , only 
l m less than HSs. We adopt the bulk heating scenario as a 
limiting case, and estimate CS evolution with the simplify- 
ing assumption of adiabatic expansion, where the radiated 
luminosity is neglected in the energy budget. The orbital 
energy, E = Eq — AE^, then depends only on R and on 
the initial orbital energy E$. The adiabatic assumption 
fails for some combinations of orbital parameters, when 
L t = L+(R 2 - \)>AE t /P. An adiabatic CS is disrupted 
earlier than a radiative one because it expands more effi- 
ciently, and because the R 5 dependence of AE t (Eq. 2) 
accelerates the heating rate. A radiative CS evolves more 
slowly, and so takes longer to reach its high luminosity 
phase, but then spends more time there. Therefore, adia- 
batic bulk heating provides a lower limit on CS lifetimes. 

The adiabatic evolution of a CS, derived from 
dAEjdt = [3(dR/dt)/R 2 = AE t /P(R), is 

t(R 2 ) - t(Ri) = Q(b,e )[f( Xl ) - f(x 2 )} , (11) 

where £q = —E /{3, x = (1 + £o)R, and where 



256V27rm 3 / 2 b 6 (l + £o ) 9 / 2 
63^/3(1 + m)T 2 (63/2) ' 



(12) 



,, , -7-10a;-16a; 2 -32.T 3 -128x 4 + 256a; 5 

fix) = . . (13) 

256y/x 9 (x - 1) 

Neglecting tidal disruption, the star expands to infinity in 

F = e(Mo)[/(l+eo)-l]. (14) 

Since the latter stages of the expansion are rapid, to es- 
timates the total time to disruption for b > 1.5. Conser- 
vative upper limits on the terminal values of the orbital 
parameters at disruption are (eo <C 1), 



mb/[20(b-l)] 
1-2/3 (6-l)/m 2 /3 

P ~ 7T {mb/[p{b - 1)]} 3/2 /y/2(l + m) 



e ~ 



(15) 



3. SQUEEZARS IN THE GALACTIC CENTER 



At a distance of 8kpc (Reid 1993), the - 3 x 10 6 M 
MBH in the GC (Ghez et al. 2000; Schodel et al. 2002) 
is the nearest and most accessible MBH. Although it is 
heavily reddened (A^-~3 m , Blum et al. 1996), deep high 
resolution IR observations provide information on the lu- 
minosity, temperature and orbits of thousands of stars near 
the MBH (Eckart et al. 1999; Figer et al. 2000; Gezari 
et al. 2002). The squeezar formation rate in the GC is 
r ~ 5 x 10 _6 yr _1 , only ~ 0.05 of the prompt tidal dis- 
ruption rate (AH03). However, squeezars are relatively 
long-lived, and so on average n~ 0.1-1 squeezars orbit the 
MBH at any given time. The leading (shortest period) 
squeezar has typically completed t/to ^n/(n+ 1) of its 
lifetime, where 7 is the mean inspiral time. 

The typical properties of the leading squeezar can be de- 
rived by averaging over all initial orbits, weighted by the 
probability of successful inspiral, which falls with increas- 
ing periapse and initial orbital period (AH03). For 1 M 
hot squeezars in the GC, this yields Lt ~ 170 Lq, T g ~ 
19000K, Aif~-2.25, P-3.6xl0 3 yr and l-e~2xl0~ 5 . 
Figure 1 shows the evolution of a HS up to disruption. 
The CS orbital bounds (Eq. 15) are similar. Note that 
the luminosity and temperature may be much higher if n 
was even modestly under-estimated. 

4. DISCUSSION 

The existence of squeezars near a MBH is a consequence 
of the flow of stars into a mass sink. We approximated the 
complex evolution of tidally heated stars by two simple 
models that likely bracket the response of real stars and 
can serve as reference for future work. 

This study of squeezars is motivated by their potential 
to probe observationally two important physical processes: 
MBH growth by tidal disruption and the effects of strong 
tides on stars. Their advantage is that they last orders of 
magnitude longer than either a tidal flare or a tidal capture 
event in stellar binaries (e.g. McMillan et al. 1987). 

It is unclear whether squeezars are spectrally distinct 
from normal stars, but their orbits should stand out since 
they require an extremely small periapse. The colors of 
stars on highly eccentric orbits can differentiate between 
surface and bulk heating, which are both plausible. Two 
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early type stars on orbits with r p ^ 120 AU, e = 0.87 and 
P=15.2 yr (Schodel et al. 2002) and r p = 60 AU, e = 0.98 
and P = 60 yr (Ghez et al. 2003) were already discovered 
in the GC. Although b ~ 35-70 is still > 10 too large for 
a squeezar, their detection shows that squeczar searches 
are feasible with deep (if ~ 18 m , cf Fig. 1) high resolution 
adaptive optics (Genzel et al. 2003). A few stars on tight 
eccentric orbits around the MBH are to be expected. For 
a simple GC model with isotropic velocities (AH03), we 
estimate that there should be ~ 120 stars on r p < 120 AU, 
P < 60 yr orbits. We note that n could be substantially 
under-estimated, since several possible scattering mecha- 
nisms were not included in the modeling of AH03. This 
possibility further motivates the search for squeezars in 
the GC. If n^> 1, a cluster of squeezars in a nearby galaxy 
may also be observable. 
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Fig. 1 — The evolution of a 1 M® hot squeezar (a — OA, Lo — 
10 4 ) in the GC, initially deflected into a 6= 1.5, P = 1.4 x 10 4 
yr orbit (to = 4.9x10 yr), as calculated analytically (full lines) 
and orbit by orbit (dashed lines). At disruption £ = 3.7x10'' yr, 
Lt = 642L , P = 210yr, l-e = 2.3xl(T 4 . The leading squeezar 
properties are given in terms of n, assuming io~to (top axis). 



REFERENCES 



Alexander, T. & Hopman, C. 2003, ApJL, in press (AH03) 

Alexander, T. & Kumar, P. 2001, ApJ, 549, 948 

Alexander, T., & Livio, M. 2001, ApJ, 560, L143 

Blum, R. D., Sellgren, K., & DePoy, D. L. 1996, ApJ, 470, 864 

Cox, J. P., & Giuli, R. T. 1968, Principles of Stellar Structure, (New 

York: Gordon & Breach) 
Eckart A., Ott, T., & Genzel, R. 1999, A&A, 352, L22 
Figcr, D. F., et al. 2000, ApJ, 553, L49 
Frank, J., & Rees, M. J. 1976, MNRAS, 176, 633 
Genzel, R., et al. 2003, ApJ, submitted 

Gczari, S., Ghez, A.M., Becklin, E.E. , Larkin, J., McLean, I.S., & 

Morris, M. 2002, ApJ, 576, 790 
Ghez, A. M., et al. 2003, ApJ, submitted 

Ghez, A. M., Morris, M., Becklin, E. E., Tanner, A. & Kremenek T. 
2000, Nature, 407, 349 



Hut, P. 1980, A&A, 92, 167 

Kumar, P., & Goodman, J. 1996, ApJ, 466, 946 
Kumar, P., & Quataert, E. J. 1998, ApJ, 493, 412 
Magorrian, J, & Tremainc, S. 1999, MNRAS, 309, 447 
McMillan, L. W., McDcrmott, P. N., & Taam, R. E. 1987, ApJ, 318, 
261 

Novikov, I. D., Petchik, C. J., & Polnarev, A. G. 1992, MNRAS, 255, 
276 

Podsiadlowski, P. 1996, MNRAS, 279, 1104 

Press, W. H., & Teukolsky, S. A. 1977, ApJ, 213, 183 

Ray, A., Kembhavi, A. K., & Antia, H. M. 1987, A&A, 184, 164 

Rees, M. J. 1988, Nature, 333, 523 

Reid, M. J. 1993, ARA&A, 31, 345 

Schodel, R., ct al. 2002, Nature, 419, 694 



